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We report on recent results on higher twist contributions to the unpolarized structure func- 
tions F£' d (x, Q 2 ) at N 3 LO in the large x region and constraints on the twist-3 contribution 
to polarized structure function g2(x,Q 2 ). 

1 Introduction 

Higher twist terms contribute to the nucleon structure functions at lower scales Q 2 . The range 
in which these terms may be safely neglected against the leading twist contributions, partly 
depends on the size of the experimental errors in the respective measurement. Highly precise 
data at low values of Q 2 allow to access these contributions, the detailed physical understanding 
of which is presently still in an early stage. It has been outlined in Refs. [2(2] how the higher 
twist contributions can be extracted in a phenomenological way in case of the structure functions 
F2(x, Q 2 ) and Fl(x, Q 2 ) in the valence quark region. In this note we report on recent results of 
an improved analysis. Another interesting question concerns the structure function g2(x,Q 2 ) 
in the polarized case, which has been measured to a higher precision during the last years [3J. 
Here we try to extract first information on the twist-3 contributions to g2(x, Q 2 )- 

2 Higher Twist Contributions to F% ,d (x, Q 2 ) 

We have carried out a QCD analysis in the valence region including more recent data from 
JLAB following earlier work pQ. In the present analysis tails from sea-quarks and the gluon 
in the valence region were dealt with based on the ABKM distributions [I]. Both the valence 
quark distributions xu v (x, Qq) and xd v (x, Qq) at Qq = 4 GeV 2 are effected only very little. 
The values of a s (M§) change marginally w.r.t. the earlier analysis P. We obtain : ct: s (M§) = 
0.1148 ± 0.0019 NLO, = 0.1134 ± 0.0020 NNLO; 0.1141 ± 0.0021 N^LO*. Here, the N 3 LO*- 
analysis accounts for the three-loop Wilson coefficients and a Pade-model for the non-singlet 
four-loop anomalous dimension, to which we attached a ±100% uncertainty, cf. p] for details. 
Furthermore, we found that the response of the individual deep-inelastic data sets in the valence 
region respond stable values, which are in accordance with the central value obtained moving 
from NLO to N 3 LO*. The present result agrees very well with determinations of a s (M§) 
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in Refs. [HE], see also [7]. A survey on the current status of a s (M§) based on precision 
measurements in different reactions has been given in [8]. In the present analysis we obtain a 
lower value of a s than the world average, cf. [8], and values being obtained in [9j[10] at NNLO. 
Reasons for the difference to the values given in [HHUI] have been discussed in Refs. [51[7] in 
detail. In particular, the partial response of a s in case of the BCDMS and SLAC data in [5UTU] 
turns out to be partly different comparing to the results in . There are also differences 
between the analyses [5] and [TU] w.r.t. several data sets contributing. 

The higher twist contributions can be determined by extrapolating the fit-results at leading 
twist for W 2 > 12.5 GeV 2 to the region 4 < W 2 < 12.5 GeV 2 ,Q 2 > 4 GeV 2 , cf. 0[n]0 The 
results for the coefficients C^ T (x) 



F 2 (x,Q 2 ) = F 2 {x,Q 2 ) 



0™[F 2 (x,Q 2 )} , Cwr(x) 



F 2 {x,Q 2 



Q 2 [GeV 2 ] 



(1) 



are shown in Figure 1, where we averaged over the respective range in Q 2 . We applied the 
target mass corrections |T3] to the leading twist contributions H The result for the higher twist 
coefficients for proton and deuteron targets depends on the order to which the leading twist 
distribution is described. The higher twist terms become smaller moving from NLO to N 3 LO*. 
Within the present theoretical and experimental accuracy the curves stabilize for x < 0.65, 
while at larger values there are still differences. 
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Figure 1: The empiric higher twist contributions to F2' d (x,Q 2 ) in the valence region, Eq. (1), extracted 
by calculating the leading twist part at NLO, NNLO, and N 3 L0*, fTT| . 



1 In Ref. [6_ also higher twist contributions for x below the valence region have been determined. 
2 An unfolding of the target mass corrections of the DIS world data for Fi and Fjj including the JLAB data, 
has been performed in 1131 recently. 
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3 g™(x,Q 2 ) 



Higher twist contributions to the polarized structure function gi(x,Q 2 ) have been studied in 
Refs. [T5J[TS] in phenomenological approaches aiming on the twist-4 contributions. However, 
the structure function g 2 (x, Q 2 ), together with other polarized electro- weak structure functions 
[TMTB] . receives also twist-3 contributions. gi{x, Q 2 ) obeys the Burkhardt-Cottingham relation 

m 

f dxg 2 {x,Q 2 ) = . (2) 
Jo 

Since the Wandzura-Wilczek relation [3D] implies, that the first moment of the twist-2 part 
vanishes separately also 

f dxg^{x : Q 2 ) = Q (3) 
Jo 

holds. The errors on the present world data from E143, E155, HERMES and NMC [3] on 
<?2(#, Q 2 ) are still large but yet one may try the fit of a profile in x. 
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Figure 2: The twist-3 contributions to g2(x,Q 2 ) subtracting the twist-2 part according to the Wandzura- 
Wilczek relation [20J using the result of [15j for the twist-2 contribution to gi(x,Q ) on experimental data 
from E143, E155, and HERMES [3] fitting the shape © (full line). Open symbols refer to data in the 
region Q 2 < 1 GeV 2 . The dashed line shows the result of a calculation at Q 2 = 1 GeV 2 given in |21| . 



In Ref. [21j the parameterization 



9T A {x) = A 



\n{x) + {1 - x) + -{I - x) 2 



{l-xf [B - (7(1 - x) + D(l 



(4) 
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has been proposed. Since the data points are measured at different values of Q 2 an evolution 
has to be performed to a common scale. Furthermore, the target mass corrections |18j have 
to be taken into account. In Figure 2 the results of the fit to g%" 3 (x,Q 2 ) are presented for 
Q 2 = 3 GeV 2 . We limited the analysis to the region Q 2 > 1 GeV 2 . The present errors are still 
large and the data of E155 dominate in the fit. We may compare with a theoretical prediction 
given in [2T]. Indeed both results are quite similar. The twist-3 contribution to the structure 
function gi(x, Q 2 ) can be obtained from that to ga(x, Q 2 ) by the integral-relation [TB] 



9? 3 



(x,Q 2 ) 



Ax 2 M 2 



9? 3 



dy 
V 



32 w3 



(5) 



cf. [22]. Due to the large errors of the data the present results are of more qualitative character. 
To study the twist-3 contributions both to the structure functions g2(x,Q 2 ) and g\{x,Q 2 ) in 
detail, a high luminosity machine, like the planned EIC [23], is needed. 



4 Conclusions 

We performed a re-analysis of the present deep-inelastic world data on proton and deuteron 
targets for the structure function F2(x, Q 2 ) in the valence region x > 0.3 accounting for remain- 
ing non- valence tails, which were calculated using the ABKM09 distributions 4 . We obtain a 
slightly lower value of a s (M^) than in our previous analysis [TJ at N 3 LO*, however, far within 
the la error range. Very stable predictions are obtained going from NLO to N 3 LO*, both for 
the valence distribution functions and a s (Af|). The values being obtained for the different sub- 
sets of experimental data in the present fit are well in accordance with our global result. We do 
not confirm the significant differences reported by MSTW between the SLAC ep and ed data at 
NNLO 0. We also disagree with the large value of NNPDF [TU] for the BCDMS data at NLO, 
which also contradicts the corresponding result by MSTW [3]. Our results are in agreement 
with those of the G JR collaboration [5, and the singlet analyses [4]|6] . We obtained an update 
of the dynamical higher twist contributions to (x, Q 2 ) in the valence region, which depends 
on the order to which the leading twist contributions were calculated. The effect stabilizes 
including corrections up to N 3 LO* in the range 0.3 < x<0.65. At larger values of x still higher 
order corrections may be needed. A first estimate on the quarkonic twist-3 contributions to the 
polarized structure function g2(x, Q 2 ) is given in a fit to the available world data on g^{x, Q 2 )- 
The contributions to gi(x,Q 2 ) are obtained by an integral relation, cf. Ref. |18) . 
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